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Abstract 

We have calculated the large-g expansion for the specific heat at the 
phase transition point in the two-dimensional g-state Potts model to the 
23rd order in using the finite lattice method. The obtained series 

allows us to give highly convergent estimates of the specific heat for q > 4 
on the first order transition point. The result confirm us the correctness 
of the conjecture by Bhattacharya el al. on the asymptotic behavior of 
the specific heat for q — > 4 + . 
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The g-state Potts model on the square lattice has a first order phase transi- 
tion for q > 4. The phase transition point and many quantities at the transition 
point are known exactly, including the free energy the internal energy |EJ and 
the correlation length^]. On the other hand, important quantities such as the 
specific heat and the susceptibility are not solved exactly. The correlation length 
increases to infinity as q — ► 4 + , so it becomes very difficult in general to eval- 
uate the quantities that are not solved exactly for q as 'small' as 5, where the 
correlation length reaches as large as several thousands. 

Here we concentrate on the specific heat at the transition point. Many meth- 
ods have been used to evaluate it including the Monte Carlo simulations || , the 
low-(and high-)temperature expansion^, and the large-g expansion^, among 
which the large-g expansion appears very promising. Bhattacharya et al. cal- 
culated the large-g series for the specific heat and higher energy cumulants to 
order 10 in z = 1/y/q- Making the Pade analysis of the series assuming the 
asymptotic form in the limit of q — > 4 + , they obtained very precise values of 
the specific heat for q > 7, whose accuracy seems to be higher than the Monte 
Carlo simulations. 

In this paper we will enlarge the large-g series for the specific heat at the 
transition point to order 23 in z. We use the finite lattice method ||, 0, |8|, 
which has been used mainly to generate the low- and high-temperature series in 
statistical systems and the strong coupling series in lattice gauge theory. The 
method can give in general longer series than those generated by the graphical 
method in lower space (and time) dimensions. In the graphical method, one 
has to list up all the relevant diagrams and count the number they appear. In 
the finite lattice method we can skip the job and reduce the main work to the 
calculation of the expansion of the partition function for a series of finite size 
lattices, which can be done using the straightforward site-by-site integration^ 
without the graphical technique. 

The model is defined by the partition function 

Z = ^exp(-/3 J ff), H = ~Y, 6 w W 

{Si} (id) 

where represents the pair of nearest neighbor sites and Sj = 1,2,- • - ,q. 

The phase transition point (3 t is given by exp (/3 t ) — 1 = ^fq. We will consider 
the free energy density in the disordered phase, which is given by 

F{fi) d = lim {L x L y )-HnZ(L x ,L y ) , (2) 

L x ,L y — >oo 

where the partition function for the L x x L y lattice should be calculated with 
the free boundary condition corresponding to the disordered phase. The large-q 
expansion of the partition function can be given through the Fortuin-Kasteleyn 
representation [^0| as 

Z(L x ,L y ) = q L ' L ^Y t N ^i^- 1 ) l 1~ i 

ho 
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= q L*L y N U3 Y l z 2j ' 1 , (3) 

1,0 

where Ni j is the number of configurations of / bonds connecting the nearest 
neighbor sites on the L x x L y lattice with L x L y — j independent clusters of sites 
and Y = [eP — l)/y/q. (Two sites connected to each other belong to the same 
cluster.) 

We define H(l Xl l y ) for each l x x l y lattice (l Xl l y — 1, 2, 3, • • •) as[|| 

H(l x ,l y ) = \n[Z(l x ,l y )/q 1 ^} , (4) 

where Z(l x ,l y ) is the partition function with the free boundary condition, and 
define W(l x ,l y ) recursively as 

W(l x ,l y ) = H(l x ,l y ) 
(') 

E (l X -l' x + my-ly + m(l' X J'y) ■ (5) 

Here the (') indicates that a term with l' x = l x and l' y = l y should be excluded 
in the summation. Then the free energy density defined by Eq.(|^) is given by 

F(p) d = ln(q) + J2W(l x ,l v ) ■ (6) 

We can prove [|| that the Taylor expansion of the W(l x , l y ) with respect to z 
and Y includes the contribution from all the clusters of polymers in the standard 
cluster expansion that can be embedded into the l x x l y lattice but cannot be 
embedded into any of its rectangular sub-lattices. Each cluster that contributes 
to the lowest order term of the W(l x , l y ) consists of a single polymer and it has 
the order of z ^+'b^ 2 . An example of such a single polymer is shown in Fig. 1. 
Therefore to obtain the series to order z N we have only to take into account all 
the rectangular lattices that satisfy l x + l y — 2 < N. If we set Y = 1 + 2/fyl, 
then we have only to keep the expansion with respect to y to order y n to obtain 
the n-th energy cumulant at the phase transition point as, 



f ln 



(We note that ^4 = (1 + y + z)j- and y = at = /3 t .) We can also calculate 



h 

duality relation 



(n) 

the series for the energy cumulant Fo at p t in the ordered phase by using the 



Ff) - F m = - g[F M - . ( 8 ) 

We have calculated the series to order N = 23 in z for n = 0, 1 and 2. The 
obtained series for the zeroth and first cumulants (i.e. the free energy and the 
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internal energy) agree with the expansion of the exactly known expressions. The 

(2) 

series for the second cumulants are listed in Table 1. The coefficients for Fo 
agree with those by Bhattacharya et al. to order 10. 
The latent heat C is known|jj to vanish at q — > 4 + as 

£ ~ 3TTX- 1 / 2 , (9) 

with x = exp(n 2 /29) and 2 cosh 9 = Jq. Bhattacharya et ai.[^2j made the 
conjecture that ' Q will diverge at q — ► 4 + as 

ifj ~ . (10) 

The constant a should be common for the ordered and disordered phases from 
Eqs. (@) and (0). Here we follow this conjecture, then the product F^C 2 can 
be expected to be a smooth function of 6, so we apply the Pade approximation 
to this quantity as 

F^C 2 = zP M (z)/Q L (z) + 0(z M+L + 2 ) , 

F^C 2 = z 2 R M (z)/S L (z) + 0(z M + L + 3 ), (11) 

where Pm{z) and Ql{z) (Rm(z) and SV^z)) are the M-th and L-th order poly- 
nomials with M + L + l < N (M + L + 2 < N). We give in Table 2 the values of 
the specific heat Cd.o evaluated from these Pade approximants for some values 
of q and present in Fig. 2 the behavior of the ratio of the F^ to x plotted versus 
9. The averages and errors are taken from all the [M, L] Pade approximants 
with M > 8 and L > 8, excluding that whose denominator has zero at some 
point in 4 < q < oo. We have checked that the duality relation (||), which is 
not respected exactly by the Pade approximants, is really satisfied within the 
accuracy for all the range of q > 4. These estimates are more precise by three 
or four orders of magnitude than (and of course consistent with) the previous 
result for q > 7 from the large-g expansion to order z w by Bhattacharya et al.^ 
and the result for q = 10, 15, 20 from the Monte Carlo simulations carefully done 
by Janke and Kappler[^|. What should be emphasized is that we obtained the 
values of the specific heat in the accuracy of about 0.1 percent at q — 5 where 
the correlation length is as large as 25000. As for the asymptotic behavior of 
F^ at q — > 4 + , the Pade data of F?^ jx and Fo /x have relatively large errors 
of a few percent arround q = 4, but their behaviors shown in Fig. 2 are enough 
to convince us that the conjecture ( p!o|) is true with 

a = 0.073 ±0.002 . (12) 

The extension of the large-g expansion to the higher energy cumulants and 
the magnetization cumulants is rather straightforward and now in progress. 
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Tabic 1: The large-q expansion coefficients i' for the second energy cumulant. 



m 


a$ (disordered) 


a™ (ordered) 











1 


2 





2 


14 


16 


3 


26 


34 


4 


118 


114 


5 


250 


254 


6 


894 


882 


7 


1936 


1944 


8 


6160 


6128 


9 


13538 


13550 


10 


39774 


39698 


11 


88360 


88360 


12 


245188 


245036 


13 


547468 


547356 


14 


1457976 


1457784 


15 


3264012 


3263316 


16 


8410284 


8410596 


17 


18868858 


18865590 


18 


47391870 


47395762 


19 


106180532 


106166828 


20 


261607968 


261629456 


21 


586199668 


586145660 


22 


1415497756 


1415594740 


23 


3174285456 


3174081000 
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Tabic 2: The specific heat for some values of q. 



q 


c d 


Co 


5 


2889(2) 


2886(3) 


6 


205.93(3) 


205.78(3) 


7 


68.738(2) 


68.513(2) 


8 


36.9335(3) 


36.6235(3) 


9 


24.58761(8) 


24.20344(7) 


10 


18.38543(2) 


17.93780(2) 


12 


12.401336(3) 


11.852175(2) 


15 


8.6540358(4) 


7.9964587(2) 


20 


6.13215967(2) 


5.36076877(1) 


30 


4.2989934145(6) 


3.4128952554(3) 
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Figure 1: An example of the cluster consisting of a single polymer that con- 
tributes to the lowest order term of the W(l x ,l y ) with l x = 4 and l y — 5. The 
closed circles are the sites, and the solid lines and the crosses are the bonds 
connecting and disconnecting the nearest neighbor sites, respectively, in the 
Fortuin-Kasteleyn representation. 
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Figure 2: The ratio of the to x plotted versus 9. The dashed and dotted 
lines represent the errors for the ordered and disordered ones, respectively. 
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